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ABSTRACT 

The  frequency  response  was  determined  for  a  second  order  servo- 
mechanism  with  backlash  in  the  system.   Frequency  response  for  sinu- 
soidal inputs  were  obtained  by  digital  computer  solutions  of  the 
equations  of  motion  of  the  system  which  had  inertia  and  viscous 
friction  on  both  sides  of  the  gear  train.   System  parameters  included 
input  amplitude,  magnitude  of  backlash,  damping  ratio  and  resilience 
of  the  gear  teeth. 

The  Control  Data  Corporation  I6c4  Computer  was  utilized  and 
Fortran  programming  was  employed. 

Results  show  the  effects  of  the  amount  of  backlash  and  the 
gear  material  resilience  on  the  output  to  input  ratios.   The  exist- 
ence of  sub-harmonic  resonance  was  found  with  large  values  of  back- 
lash . 
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L.   BACKGROUND  AND  INTRODUCTION 

Most  servomechanisms  have  a  mechanical  connection  between  the 
driving  components  and  the  driven  or  positioning  unit.   Usually  this 
connection  is  a  gear  train  which,  in  most  practical  cases,  has  a 
finite  amount  of  backlash  in  the  gearing  mechanism.   This  inherently 
renders  the  servo  system  non- linear  and  requires  analysis  techniques 
which  are  based  on  approximations  o:  which  necessitate  use  of  laborious 
graphical  methods  which  may  rely  on  simplifying  assumptions.   The 
backlash  between  the  gear  teeth  cannot  always  be  eliminated  without 
possibly  introducing  other,  perhaps  less  desirable  effects  such  as 
coulomb  friction  and/or  increased  viscous  friction.  Manufacturing 
costs  may  also  prohibit  consideration  of  gears  with  no  backlash. 

The  relative  distributions  of  viscous  friction  and  inertia  on 
each  side  of  the  gearing  have  direct  effects  on  the  behavior  of  the 
load.   Graphical  phase  plane  and  analog  simulation  analysis  techniques 
do  not  readily  adapt  to  practical  consideration  and  adjustment  of 
these  distributions.  Also,  they  cannot  handle  the  case  of  non- 
plastic  impact  of  the  gear  teeth.   This  has  led  to  the  solution  of 
this  problem  by  the  use  of  digital  computers. 

Several  investigations  have  been  made  of  second  order  servo 
behavior  using  the  digital  computer  at  the  U.  S.  Naval  Postgraduate 
School.   These  include,  limit  cycle  investigation  by  New  for  plastic 
impact  (Ref.  \)   and  by  Anderson  and  Luckett  for  plastic  and  elastic 
impact  (Ref.  2),  and  a  transient  response  investigation  for  plastic 
impact  by  Andrews  and  Kelley  (Ref.  3)- 
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This  thesis  extends  the  coverage  of  this  problem  to  a  steady- 
state  frequency  response  investigation,  with  both  plastic  and  elastic 
impact,  for  a  periodic  (sinusoidal)  input  of  varying  frequency  and 
amplitude.   The  general  approach  to  the  investigation  was  similar  to 
that  used  by  Anderson  and  Luckett.   The  equations  of  motion  of  the 
system  were  normalized  so  behavior  would  be  independent  of  the  un- 
damped natural  frequency  of  the  system.   The  parameters  inertia, 
viscous  friction,  elastic  bounce   coefficient  and  damping  ratio 
were  in  the  same  ranges  of  values  used  by  Anderson  and  Luckett  so 
that  a  correlation  of  results  could  be  made  if  desired. 

Computations  were  made  on  a  CDC  1604  digital  computer  using  a 
Fortran  System  program  which  had  definite  advantages  over  the  machine 
language  programs  used  in  the  previous  investigations  mentioned  above. 
Solution  of  the  system  equations  was  based  on  the  phase  plane  approach 
as  to  whether  the  system  was  operating  combined  or  separated.   Even 
though  the  system  was  normalized,  the  computer  program  was  made  gen- 
eral to  permit  extension  of  the  values  of  the  parameters  if  desirable. 

This  thesis  includes  investigation  of  the  effects  of  the  amount 
of  backlash,  the  amplitude  of  the  periodic  input  and  the  amount  of 
system  damping  for  both  plastic  and  elastic  impact  of  the  gear  teeth. 
The  only  non-linearity  considered  was  backlash. 

The  total  system  was  piecewise  linear  in  that  one  linear  differen- 
tial equation  described  the  motion  of  the  system  when  the  driving  motor 
was  combined  with  the  load  and  two  different  linear  differential  eq- 
uations described  the  separated  motions  of  the  motor  and  load. 
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Boundary  conditions  for  the  differential  equations  were  based  on  phase 
plane  slopes  and  on  equations  based  on  the  laws  of  conservation  of 
energy  and  momentum  which  were  related  to  whether  gear  contact  was 
plastic  or  elastic. 

The  original  intention  was  to  determine  the  effects  of  different 
distributions  of  inertia  and  viscous  friction  between  the  driving  unit 
(motor)  and  the  driven  unit  (load).  However,  because  of  the  great 
amount  of  digital  computer  operating  time  required,  it  was  necessary 
to  confine  the  investigation  to  one  ratio  of  inertias  and  one  ratio  of 
frictions . 

Results  from  the  digital  computer  calculations  provided  data  from 
which  the  magnitude  of  the  output  and  the  phase  relation  between  input 
and  output  could  both  be  determined.   This  thesis  is  concerned  only 
with  the  magnitude  of  the  output  as  compared  with  amplitude  of  the  in- 
put sinusoid. 

Special  effects  or  phenomena  such  as  jump  resonance  and  sub- 
harmonic  resonance,  which  are  known  to  exist  under  certain  conditions, 
were  looked  for  in  the  course  of  the  investigation. 

In  the  following  sections  the  system  equations  are  developed  and 
the  evolution  of  the  computer  program  is  described,  A  print  of  one 
version  of  the  Fortran  language  computer  program  and  the  computer 
solution  flow  chart  are  included  as  appendices . 


2.   SYSTEM  DESCRIPTION  AND  DEVELOPMENT  OF  EQUATIONS 


2  .1  Description  of  System 

For  the  purpose  of  this  investigation,  a  second  order  system  was 
chosen.   The  block  diagram  of  the  system  is  as  follows: 
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Fig.   2-1 


This  system  is  the  same  system  that  was  reported  by  Anderson  and 
Luckett  in  Ref.  2.   Except  for  the  non-linearity,  the  system  is  a 
simple  second  order  servomechanism.   There  is  only  one  non-linearity 
present,  and  that  is  backlash.   The  backlash  is  present  in  the  gear 
train  which  separates  the  motor  from  the  load.   The  effect  of  the 
backlash  is  to  allow  movement  of  the  load  and/or  the  motor  inde- 
pendent of  the  movement  of  the  other,  so  long  as  this  movement  does 
not  exceed  the  amount  of  backlash.  Also  included  as  part  of  the  gear 
train  is  the  gear  ratio;  however,  this  has  no  special  effect  on  the 
system  other  than  is  normally  expected. 

Separated  by  the  backlash  is  the  motor  inertia  and  friction  and 
the  load  inertia  and  friction.  The  distribution  of  these  quantities 
cannot  be  changed  in  the  describing  equations  of  the  system  due  to 


the  backlash.   In  the  linear  case,  the  inertia  and  friction  can  be 
lumped  on  either  side  of  the  gear  train  without  affecting  the  system 
operation,  but  this  cannot  be  done  when  there  is  backlash  present. 
The  motion  of  the  system  can  be  described  by  a  combination  of 
three  linear  differential  equations.  When  the  gear  teeth  are  touching, 
a  single  differential  equation  describes  the  motion .   Once  the  gear 
teeth  separate,  one  linear  differential  equation  describes  the  motion 
of  the  motor,  and  another  describes  the  motion  of  the  load.   Upon 
recombining  of  the  gear  teeth,  the  single  differential  equation  again 
describes  the  motion.  At  the  instant  of  separation  or  recombination, 
the  appropriate  differential  equations  have  as  initial  conditions  the 
last  values  of  the  equation  which  described  the  motion  prior  to  the 
separation  or  recombination.   The  change  over  from  one  describing 
equation  to  the  pair  of  equations  and  vice  versa  will  be  discussed 
further  in  the  equation  development  section. 

2.2  Assumptions 

Prior  to  developing  the  equations  which  describe  the  motion  of  the 
system  under  consideration,  the  assumptions  which  have  been  made  will 
be  presented.   These  assumptions  are  essentially  those  made  in  Ref .  2 
and  are  repeated  for  convenience  but  with  some  modification. 

The  following  assumptions  were  made: 

1.  The  gear  teeth  were  initially  in  contact,  and  the  initial  con- 
ditions of  the  system  were  all  equal  to  zero. 

2.  Plastic  deformation  of  the  gear  teeth  during  steady  state  con- 


tact  or  on  impact  and  any  torsional  deformations  of  driving 
shafts  are  negligible. 

3.  The  inertias  of  the  gears  and  drive  shafts  are  considered  as 
part  of  the  load  or  motor  inertia  depending  on  their  attach- 
ment in  the  system. 

h .      The  law  of  conservation  of  energy  is  completely  satisfied  in 
the  perfect  elastic  case,  and  the  law  of  conservation  of 
momentum  is  satisfied  in  the  imperfect  elastic  cases  includ- 
ing the  plastic  case.   When  only  the  law  of  conservation  of 
momentum  is  required  to  be  satisfied,  the  energy  lost  by  the 
system  is  dissipated  in  the  form  of  heat  by  infinitesimal 
deformations  of  the  gear  teeth, 

5.  The  gear  teeth  are  in  contact  only  instantaneously  during  im- 
pact for  all  cases  except  the  plastic  case,  and  the  impulse 
torques  of  drive,  friction,  bearing  supports,  etc.,  are  zero 
during  impact . 

6.  The  coefficient  of  restitution  of  the  two  opposing  gear  teeth 
is  the  same  or  can  be  described  by  an  equivalent  coefficient 
if  the  gear  teeth  are  of  unequal  coefficients. 

7.  Backlash  is  assumed  to  be  equal  at  all  points  on  the  gear  cir- 
cumference, and  the  backlash  is  measured  at  the  output  shaft. 

2 .3   Equation  Development 

From  Fig.  2-1,  the  error  signal,  E,  is  the  driving  signal.   This 
is  equal  to: 


E  =  eR  -ec  (2-D 


The  torque  developed  by  the  motor,  T  ,  is  equal  to: 


Tm  =ke  -k  (eR-ec)  ^s) 


Equations  (2-1)  and  (2-2)  are  true  when  the  load  and  motor  are  combined 
and  when  they  are  separated. 
The  load  torque,  T  ,  is: 

Id 

tl  =  JLec+fLec  <*3> 

When  the  load  and  motor  are  combined,  the  load  torque  when  referred  to 
the  motor  side  of  the  gear  train  is: 

tlm=/>tl  =/>(jLec  +  fLec)  (2--o 

The  forward  or  positive  direction  of  the  system  is  defined  to  be 
when: 

ec  =/>  eM  (2-5) 

and  the  backward  or  negative  direction  when: 

©C  "/^M  +  A  (2-6) 

By  differentiating  either  equation  (2-5)  or  (2-6)  with  respect  to 
time  yields  the  same  results.   These  are: 


ec  v  Qm  (2-?) 

and 

Substituting  equations,  (2-7)  and  (2-8),  into  equation  (2-U) 
yields  after  rearranging: 

T     -p2fj  9  +f  fl)  (2-9) 

1  LM    /    \JL  WM  '  I  L  ^tAj 
The  driving  torque  of  the  motor  must  equal  the  load  torque  referred  to 
the  motor  side  (combined  system  only)  and  the  load  torque  of  the  motor 
itself. 

T      -     j    C\     +  f '    f)     +  T  (2-10) 

Substitution  of  equations  (2-1),  (2-7),  (2-8),  and  (2-9)  into  (2-10) 
after  simplifying  gives: 


de  +(JM  V>*JJ  0C  +  (J.  vU)  9<  "  (JM  v»JL)  9r     (     } 


Rewriting  equation  (2-11) 

0c  -h  2  J0Jn  Qc+   60*  0C   =  60*    0R  (2-12) 


where 
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z         K/> 

<^n  =   ,  /  z    |  (2-13) 


and 


2^^  =  -r^T-  (S-1U) 

Equation  (2-12),  therefore,  defines  the  motion  of  the  system  when  the 
gear  teeth  are  combined. 

Whenever  the  gear  teeth  are  separated,  there  is  no  driving  force 
on  the  load,  and  its  motion  is  defined  by  the  following  equation: 

©L  +-^0L  =  0  (2-15) 

To  be  noted  at  this  point  is  the  change  in  notation  for  the  load. 
In  order  to  distinguish  between  combined  and  separated  operation,  the 
subscript  C  or  L  is  used  respectively  for  the  output. 

Combined  operation: 


Output  position 

»c 

Output  velocity 

Qc 

Output  acceleration  - 

a 

Separated  operation: 

Output  position 

eL 

Output  velocity 

a 

Output  acceleration  - 

eL 

Although  this  convention  in  notation  was  adopted,  it  should  be  remembered 
that  the  same  quantity  is  described  by  two  different  symbols,  depending 
on  whether  the  system  is  combined  or  not . 

By  the  gear  teeth  being  separated,  the  load  torque  is  not  reflected 
back  to  the  motor  side  of  the  gear  train.   Therefore,  equation  (2-10) 
reduces  to: 

J       =   I   H  j-   f  Q  (2-16) 

or  by  substituting  equation  (2-2)  into  (2-16)  and  after  rearranging 
reduces  to: 

eM  +  %  oM  +  4-  ec  =  -f  eR  (2.17) 

WM  WNA  J/v\ 

Since  this  investigation  is  directed  towards  the  frequency  response 

of  the  system,  the  forcing  function  is: 

QR  ~  A  sin  cot  (2-18) 

If  all  angular  frequencies  are  divided  by  the  undamped  natural  angular 
frequency,  od  ,  this  effectively  is  the  same  as  setting  a)  equal  to  1. 
By  this  normalization,  the  angular  frequency  in  equation  (2-18)  is 
really  u)/a)  .   Therefore,  the  results  are  applicable  to  any  simple  second 
order  system  having  any  natural  frequency  a)  ,  if  the  forcing  function 
angular  frequency  is  multiplied  by  the  cu  of  that  system. 

Equation  (2-12)  describes  the  motion  of  the  system  when  the  gear 
teeth  are  combined.   Equations  (2-15)  and  (2-17)  describe  the  motion  of 
the  load  and  motor  respectively  when  the  gear  teeth  are  separated. 
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Left  to  be  determined  now  is  the  instant  when  separation  or  combination 
occurs  . 

Use  of  phase  plane  analysis  is  necessary  to  determine  the  instant 
of  separation  when  the  gear  teeth  are  combined.  By  using  velocity  and 
position  coordinates  for  the  phase  plane,  the  slope  of  a  trajectory  is: 


and  specifically  for  the  combined  system: 


Nd  G 
=  — —  (2-19) 


M  =  44  =  -&-  (2-20) 


'c 


6  ec      0, 

With  proper  manipulation  of  equation  (2-12),  the  following  obtains: 

U       -        CUn    8r  -C0n6c  -  2  JoJn  9C  (2_21) 

©C 

Using  similar  techniques  to  find  the  load  alone  slope  on  the  phase 
trajectory  from  equation  (2-15)  yields: 

NL  -  ~4-  (2-22) 

Assuming  that  the  system  is  combined,  it  will  stay  combined  as 
long  as  the  magnitude  of  the  load  alone  trajectory  slope  is  greater  than 
the  magnitude  of  the  combined  trajectory  slope.   Symbolically,  the 
system  will  remain  combined  so  long  as  the  following  is  true: 
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NL|  >  |NC  (2-23) 

Separation  occurs  when  the  combined  trajectory  slope  equals  the 
load  alone  trajectory.  This,  therefore,  is  the  criterion  to  be  satis- 
fied for  transferring  the  motion  description  of  the  system  from 
equation  (2-12)  to  equations  (2-15)  and  (2-17). 

For  the  system  to  recombine  after  separation,  the  gear  teeth  must 
come  into  contact  with  such  velocities  that  they  do  not  bounce  apart  in 
the  elastic  case.   To  recombine  in  the  plastic  case,  the  gear  teeth 
merely  must  come  into  contact. 

In  the  general  case,  two  criterion  must  be  satisfied  for  the  load 
to  recombine  with  the  motor.  These  are: 

1.  The  gear  teeth  must  touch. 

2.  The  velocities  of  the  load  and  motor  must  be  such  that  they 
are  equal  after  impact.   (Elastic  case  only) 

Unless  equation  (2-5)  or  (2-6)  is  true,  the  gear  teeth  cannot  be 
in  contact.   Therefore,  when  either  equation  (2-5)  or  (2-6)  is  satis- 
fied, the  gear  teeth  have  come  into  contact.   In  making  this  test, 
however,  it  must  be  remembered  that  t?L  is  actually  used  in  place  of 

sOc    in  equations  (2-5)  and  (2-6).   It  is  re-emphasized  that  these 
symbols  actually  represent  the  same  quantity  but  are  obtained  by  solving 
different  equations. 

If  the  impact  is  plastic,  the  gear  teeth  will  stay  together.  How- 
ever, in  the  elastic  case,  they  do  not  necessarily  have  to  remain  in 
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contact.   During  impact,  the  law  of  conservation  of  momentum  must  be 
satisfied;  and  after  impact,  the  velocity  ot  the  load  and  motor  must 
be  the  same  for  the  teeth  to  remain  in  contact.   By  using  the  lav  of 
conservation  of  momentum,  the  respective  velocities  can  be  determined. 
In  Ref.  2,  the  coefficient  of  restitution  was  defined.   This  is: 


e  -  - 


/°6M  -9L 


(2-21+) 


or  rearranging: 


9L  -e{p%-QL 


(2-2UA) 


where  the  primed  terms  represent  the  respective  velocities  after  im- 
pact.  To  satisfy  the  conservation  of  mementum  law: 


J 


P 


f(eM-eMj=(e;-eL 


(2-25) 


must  be  true.   To  be  noted,  equation  (2-25)  ls  independent  of  the 
coefficient  of  restitution. 

Solving  equation  (2-24)  for   L7M,  gives  equation  (2-2^A);  and  sub- 
stituting the  result  into  equation  (2-25),  the  following  expression  for 


eL 


obtains 


eL  = 


J 


M 


(JL  +/>U 


pQM(e  +  \)+Q 


■SJL 


d 


-  e 


M 


(2-26) 


U 


If  the  gear  teeth  have  come  into  contact,  equations  (2-26)  and 

•  /      •  / 
(2-2^A)  can  be  solved  for  wL  and  W^  respectively.   Unless 
•  /    •  / 
W|_  =  (3)^  ,  the  gear  teeth  will  not  remain  together. 

In  summary,  assuming  that  the  system  starts  in  the  positive  direc- 
tion from  zero  with  the  gear  teeth  in  contact,  the  operation  is  as 
follows: 

1.  Equation  (2-12)  defines  the  motion  of  the  system  until  equation 
(2-23)  is  no  longer  true. 

2.  Once  equation  (2-23)  is  no  longer  true,  equations  (2-15)  and 
(2-17)  describe  the  motion,  taking  as  their  initial  values 
the  terminal  values  of  (2-12).   In  general,  the  direction  of 
the  system  must  be  determined  so  that  the  proper  value  for 

Oj/  can  be  used  to  initialize  equation  (2-17).   If  t/c  is 
positive,  equation  (2-5)  applies.   If  wc  is  negative, 
equation  (2-6)  applies.   In  transferring  from  equation  (2-12) 
to  (2-15),   0L  =  6C  . 

3.  During  separated  operation,  recontact  is  determined  from  equa- 
tion (2-5)  or  (2-6). 

h .      In  the  event  of  recontact,  equations  (2-26)  and  (2-2UA)  are 
used  to  determine  if  bounce  occurs.   (Required  only  in  elas- 
tic case) 

5.  Finally,  when  there  is  no  bounce,  description  of  the  system's 
motion  reverts  to  equation  (2-11).   This  time,  the  terminal 
value  of  t>i_  is  used  to  initialize  equation  (2-12)  along 
with  0L  . 

6.  The  steps  of  operation  continue  this  pattern. 
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2 .h      Parameters 

By  referring  to  equations  (2-12),  (2-15),  (2-17),  and  (2-24),  the 
following  parameters  are  available  for  varying: 

i.  y 

2.  LOr\ 

3.  ©R 

*.    Jm 

5.  TtA 

6.  JL 

8.      e 

By  setting  6Jn  equal  to  one,  the  system  was  made  independent  of 
its  natural  frequency.   This  normalization  was  discussed  previously. 
Not  shown  in  the  preceeding  list  as  a  parameter  is  the  gear  ratio, 
yO  ;  of  the  gear  train.   Since  the  gear  ratio  was  considered  to  have 
no  special  effect  on  the  system,  its  value  was  set  equal  to  one.  Also 
not  shown  is  the  magnitude  of  backlash.   This  was  not  used  as  a  para- 
meter; however,  it  was  assigned  a  value  other  than  unity. 

The  motor  inertia  and  load  inertia,  J_,  and  JT  respectively, 

M      L 

were  combined  such  that  with  P    equal  to  one  their  sum  was  one. 
(i.e.,  J  +  J°   2  ^T  -  1)   This  means  that  any  combination  of  J  and 
J   can  be  used,  so  long  as  their  sum  is  one. 

By  solving  equation  (2-13)  for  the  system  gain  constant,  K,  the 
following  is  obtained. 
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Since  a)  is  one,  (j  4-  y°  2  J  )  is  one  and  £>    is  one;  the  gain  con- 
stant, K,  must  also  be  one.   This,  therefore,  was  the  value  assigned  to 
K  in  this  investigation, 

Anothe-r  variable  was  defined  by  combining  the  load  and  motor 
friction,  f  and  f  ,  similar  to  the  way  the  inertias  were  combined. 
This  new  variable  called  total  friction  is  defined  as  follows: 

fT    r  fM  +/?k  (£-28) 

The  friction  ratio  of  the  system  is  defined  to  be  f  /f  . 

L-    JL 

Substituting  the  value  of  one  for  o>  ,  P   and  (j  +  /D    2  J  )  in 
equation  (2-lU),  the  following  results 


or 


J  _   f^  +  ^L  (2.29A) 

2 

Assigning  a  value  to  the  damping  factor,  *J  ,  fixes  the  sum  of  ( f  +  f  ) 

To  determine  the  value  of  fw  or  f . ,  a  choice  also  must  be  made  for  the 

M     L 

friction  ratio.   Choice  of  jf  and  the  friction  ratio,  therefore,  fixes 

f„  and  f  . 
M       L 

For  this  investigation,  the  imput,  \Jo    >  was  a  sinusoid  defined  by 
equation  (2-18).   The  angular  frequency  of  oscillation,  od,  was  effec- 
tively normalized  by  making  a)  equal  to  one.   However,  the  values 
assigned  to  CD  were  varied  over  the  desired  range  of  the  frequency 
response.   There  is  also  no  restriction  on  the  choice  of  the  forcing 
function  amplitude  A. 
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Equation  (2-24)  defines  the  coefficient  of  restitution  of  the  gear 
teeth.   This  definition  imposes  no  restriction  on  its  value;  practi- 
cally however,  it  must  be  less  than  one  and  greater  than  zero. 

Having  so  many  parameters  available,  it  was  necessary  to  restrict 
their  variation.   Listed  below  are  the  values  assigned  to  the  various 
parameters.   It  was  felt  that  this  choice  would  produce  the  most  in- 
teresting results.   The  method  used  to  vary  parameters  will  be  covered 
in  Section  3. 
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PARAMETER 
Natural  frequency 
Gear  ratio 
Backlash 
Load  inertia 
Motor  inertia 
System  gain  constant 
Motor  friction 
Load  friction 
Friction  ratio 
Damping  ratio 
Forcing  function  frequency 
Forcing  function  amplitude 
Coefficient  of  restitution 


TABLE    II- 

•  i 

Values  Used 

for 

Parameters 

SYMBOL 

VALUES   USED 

U)n 

1.0 

J° 

1.0 

•■ 

A 

.1,    .05 

Jl 

.5 

J\A 

.5 

K 

1.0 

f  M 

.12, 

.21+,    .36,    .48,    .60,    .72, 

.96 

fL 

.08, 

.16,    .24,    .32,    ,40,    .48, 

.64 

H 

.4 

3 

.05, 

.1,    .2,    .3,    .4,    .5,    .6, 

.8 

00 

.1    to    10.0 

A 

10,    1.0,    .5,    .095,    .12 

e 

0,    .6,    .8 
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3.   DIGITAL  COMPUTER  PROGRAMMING  AND  METHOD  OF  SOLUTION 

3 .1   Digital  Computer  Program 

All  of  the  frequency  response  results  were  obtained  by  com- 
putations using  the  Control  Data  Corporation  160^-  high-speed  digital 
computer.   The  computer  programming  was  done  using  the  Fortran  system 
with  program  input  either  by  magnetic  tape  or  by  punched  cards.   Com- 
puter output  was  put  onto  magnetic  tape  and  results  were  printed 
from  the  magnetic  tape  by  use  of  an  IBM  1^01  high-speed  line  printer 
unit . 

The  Fortran  system  was  used  because  Fortran  programming  permitted 
inserting  the  equations  of  the  servo  system  into  the  program  in  forms 
very  close  to  their  actual  mathematical  forms.   The  use  of  Fortran 
also  allowed  a  great  deal  of  flexibility  in  easily  modifying  the  program, 
in  changing  system  parameters,  in  inserting  input  data  and  in  extract- 
ing the  output  results  in  a  properly  labeled  decimal  format. 

A  print  of  a  typical  Fortran  source  program  used  in  this  inves- 
tigation is  shown  in  Appendix  B.   This  Fortran  source  program,  includ- 
ing sub- routines,  represents  about  1600  machine  language  instructions. 

The  system  equations  to  be  solved  by  the  computer  were  developed 
in  Section  2.   The  three  differential  equations  describing  the  motion 
of  the  combined  system  and  the  motion  of  the  separated  motor  and  load 
were  solved  by  subroutines  which  were  parts  of  the  overall  Fortran 
program.   These  sub- routines  were  based  on  the  Runge-Kutta  numerical 
method  for  the  solution  of  ordinary  differential  equations  using  a 
four  stage  process.   Each  of  the  Runge-Kutta  sub- routines  utilized  its 
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own  integral  sub-program. 

Even  though  the#  system  being  investigated  was  normalized  (a)  =  1) 
and  a  unity  gear  ratio  was  chosen  {yO=    I),    the  program  was  generalized 
to  provide  for  changing  03  and  yO    to  any  desired  magnitudes.   Similarly, 
it  was  possible  to  select  any  value  of  system  gain  K,  however  a  value  of 
K  =  1  was  used  for  all  calculation  for  reasons  discussed  in  Section  2. 

As  shown  in  Section  2,  the  forcing  function  or  input  to  the  system 
in  this  frequency  response  investigation  was  C7R  =  A  sin  cot.   The 
amplitude  A  of  the  sinusoidal  input  could  readily  be  set  prior  to  any 
computer  run.   The  sinusoidal  input  was  obtained  from  the  Fortran 
Reference  Library  sine  function  which  is  an  integral  part  of  the  Fortran 
system.   The  input  frequencies  cjd  for  each  frequency  response  run  were 
put  into  the  program  in  an  array  of  data.   Each  a)  was  consecutively 
inserted  into  the  calculations  by  the  program  until  all  desired  a)' s 
were  completed.  A  knowledge  of  Fortran  programming  is  necessary  to 
thoroughly  understand  how  this  is  accomplished. 

It  was  necessary  to  choose  integration  intervals  which  would 
yield  accurate  results  for  the  Runge-Kutta  solutions  of  the  differen- 
tial equations  yet  which  would  not  require  more  calculations  than 
necessary.   It  was  found  that  for  values  of  J  ^-0.1  the  following 
relationship  was  satisfactory  for  determining  the  intervals. 
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i        o.ooi  ,    ., 

Integration  interval      /\  L  —  — ; IJ--U 

CO 

This  resulted  in  2000  x  2TT  solutions  per  input  cycle  for  all  values 
of  input  frequency  o>.   The  values  of   At  were  also  put  into  the  pro- 
gram in  an  array  of  data  and  each  value  was  inserted  into  the  calcu- 
lation along  with  its  corresponding  a). 

With  combination  of  low  a)  and  j*  <C0.1,  it  was  necessary  to  re- 
duce At  to  a  value  less  than  obtained  from  equation  (3-1)  in  order  to 
maintain  sufficient  accuracy. 

Since  the  damping  ratio  J°  for  any  run  was  dependent  on  the  rela- 
tive magnitudes  of  f  and  f  ,  these  parameters  (j",  fM  and  f  )  were 
set  prior  to  each  separate  run  and  were  constant  for  that  run.   All 
runs  were  made  with  a  constant  ratio  of 

-  0.4 


Even  though  the  distribution  of  inertia  between  motor  and  load 
was  even  for  all  calculations,  the  program  provided  for  changing  this 
distribution  to  80$>  of  the  total  inertia  to  either  the  motor  or  load 
side.   This  was  built  into  the  program  and  arrays  of  values  for  J  and 
J  were  put  in  as  data  similar  to  the  arrays  of  a)  and  At. 

Li 

The  amount  of  backlash  between  the  gear  teeth  was  fixed  before 
each  computer  run;  in  fact  the  value  used  for  all  but  one  rar,  was 
A  =  0.1  radian.  A  value  of  A  =  0.05  was  used  in  a  special  run  to 
check  the  effect  of  changing  A  and  A  proportionally. 
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Because  of  the  large  number  of  solutions  per  input  cycle  (  OOCTf) 
and  since  it  was  not  necessary  to  record  all  of  these  results,  the 
program  was  designed  to  print   the  results  for  every  one  hundredth 
position  computation.   This  provided  sufficient  results  to  obtain  the 
frequency  response  of  the  system. 

Only  the  steady  state  response,  after  the  system  transients  had 
effectively  died  out,  was  desired.   Therefore  a  delay  of  five  time  con- 
stants (based  on  a  similar  linear  system)  was  included  in  the  program  so 
outputs  would  not  be   printed  until  this  transient  period  had  passed. 
The  program  mnenomic  for  this  was  DELAY. 

It  was  necessary  to  tell  the  program  how  long  to  compute  before 
switching  to  the  next  value  of  a).   This  was  done  by  causing  the  compu- 
tations to  shift  to  the  next  cu  when  a  predetermined  amount  of  problem 
time  had  elapsed  after  the  transient  period.  This  time  was  based  on  a 
desired  number  of  input  cycles.   The  program  mnemonic  for  this  was 
designated  CYC. 

3.2  Sequence  of  Solution  of  System  Equations 

The  starting  condition  for  each  frequency  response  run  was  with 
the  gear  teeth  in  contact,  i.e.  the  motor  and  load  were  combined,  with 
this  equation  describing  the  system  motion: 

0C  -f  2  Jcon  8C  +  cu^  0C  =  03n   A  5 in  cot  (3-2) 

This  differential  equation  was  solved  as  each  increment  of  sinusoidal 
input  was  added.   The  corresponding  value  of  0cwas  calculated  and  the 
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combined  phase  plane  slope  was  computed  by 


Nc  =  -^-  (3-3) 


This  slope  was  compared  with  the  computed  value  of  the  load-alone 
slope,  which  was  constant  for  the  run; 


K- 


(3-"0 


Ju 


When  H  =  N.  the  motor  and  the  load  were  commencing  to  separate.  If 
the  system  velocity  0cwas  positive  at  separation  the  gear  teeth  were 
separating  from  the  original  contact  position, 

If  0  was  negative,  the  teeth  were  separating  from  the  opposite  side  of 
the  backlash. 

After  separation,  the  equation  describing  the  motor  motion  was 

eM+-^eM+TeL=Jf  A  sin  oo I  (>5) 

WM  \J(A  ^bA 

This  differential  equation  was  solved  for  each  input  increment.   At 
the  same  problem  time  t,  the  following  differential  equation  for  the 
load  alone  was  solved: 

9L  +  ^r6L  =  0  (>6) 

From  these  solutions,  the  positions  of  the  motor  and  load  were 
compared  in  order  to  determine  when  the  gear  teeth  had  re-touched  on 
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either  side  of  the  backlash,   A  condition  to  be  satisfied  for  re- 


touching was 


°r 


eL  =/  e 


r-1 


0L=/>9+A 


(3-7) 


(3-8) 


When  the  teeth  did  make  contact,  the  load  velocity  after  touching  was 
calculated  from 


e.  = 


j 


M 


L 


JlA+Pdl 


e 


and  the  motor  was 


•  /  _  eL  -  e  (/>  eM  -  eu 
■^ p 


(3-9) 


(3-10) 


The  following  equation  was  then  solved  to  see  if  the  teeth  stayed 
combined 


is    />eM  -  eL  ^  €  ? 


(3-U) 


Where  6  was  an  arbitrarily  small  number  which  was  used  to  determine 
when  the  velocities  were  almost  equal  in  the  case  of  elastic  impact. 

When  the  teeth  stayed  combined,  the  program  returned  to  the 
combined  differential  equation  (3-2)  and  the  complete  sequence  was 
repeated . 
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3  «3   Printing  of  Results  . 

Solution  of  any  of  the  three  differential  equations  of  the  system 
produced  the  values  of  the  respective  position  G  and  velocity  © . 
Since  the  position  and  velocity  of  the  combined  system  were  defined  as 
the  load  or  output  position  and  velocity,  these  values  of  0cand0cwere 
printed  as  the  outputs  when  the  gear  teeth  were  in  contact.  When  the 
motor  and  load  were  separated,  the  load  alone  was  printed  as  the 
output,  and  the  motor  position  and  velocity  were  also  printed  for  the 
same  problem  time  t. 

For  each  print  of  system  position  and  velocity,  the  problem  time 
and  the  input  position  were  also  printed.   Prints  were  made  for  each 
one  hundredth  solution  of  the  differential  equation  or  equations  con- 
cerned.  Special  prints  of  problem  time,  position(s),  velocity(s)  and 
input  were  made  just  after  the  gear  teeth  separated  or  combined. 
A  typical  print  of  computer  output  is  shown  in  Appendix  C.   This 
shows  that  output  position  and  velocity  (©and  0  )  only  were  printed 
when  the  system  was  combined,  while  motor  and  load  positions  and 
velocities  also  were  printed  when  the  teeth  were  separated. 

Only  five  decimal  places  of  results  were  printed  because  of 
limitations  in  the  accuracy  of  plotting  results.   The  computational 
accuracy  was  somewhat  greater  since  the  Runge-Rutta  method  was 
accurate  to  about  (At)  ,  which  was  10    for  the  largest  At  used. 

A  complete  computer  program  flow  chart  and  explanation  is  in- 
cluded in  Appendix  A. 
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3 ,h      Computation  Time 

y  

Computer  operating  time  to  obtain  a  steady  state  output  varied 
considerably  for  different  values  of  ^   and  for  the  different  input 
frequencies  CO.   For  example,  to  compute  through  a  five  time  constant 
transient  period,  with  a  low  damping  ratio  of  ^f   =  0,1  (long  time  con- 
stant) and  with  an  a)  of  2.0  (small  integration  interval),  required 
about  28  minutes  of  computation.   On  top  of  this,  it  took  about  one 
and  three-quarter  minutes  to  compute  and  print  the  outputs  for  each 
input  cycle.   From  one  to  eight  cycles  of  input  (after  the  estimated 
transient  period)  were  required  to  obtain  a  steady  state  peak  output. 
For  the  other  extreme  of  a  high  J  of  0.8  and  a  low  cc  of  0.1,  the 
total  computation  and  print  time  required  to  obtain  an  output  peak 
took  only  about  two  minutes. 

It  is  apparent  that  the  investigation  required  a  great  amount 
of  computer  time.  All  computer  runs,  including  test  runs,  took  a 
total  of  about  83  hours  of  computer  operation. 

An  attempt  was  made  to  develop  a  computer  program  which  would 
approximate  the  energy  stored  in  the  system  for  steady  state  at  the 
first  input  frequency  and  would  retain  this   stored  energy  when  a) 
was  changed.   It  was  hoped  this  would  materially  reduce  the  transient 
period  in  reaching  steady  state  at  the  new  0).   Test  runs  were  promising 
but  time  did  not  permit  refinement  of  this  program  to  where  it  could 
be  used  in  the  investigation. 
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k.      RESULTS  AND  DISCUSSION 


k .  1   Determining  M  from  the  Computer  Output 

From  standard  frequency  response  nomenclature,  the  symbol  M  is 
defined  to  be  the  magnitude  of  the  ratio  of  system  output  to  system 
input.   Symbolically  this  is: 


M  = 


a 


e, 


(*.i) 


where 


0, 


=  A 


Appendix  C  shows  a  sample  print  out  of  data  from  the  computer. 


a 


One  of  the  quantities  printed  is  output  position.   To  determine 
it  was  necessary  to  pick  off  the  value  of  output  position  at  its  maximum 
or  minimum.   In  general,  there  were  two  peaks.   One  was  a  positive  peak 
and  the  other  a  negative  peak.   These  two  peaks  came  about  because  the 
input  was  permitted  to  go  through  only  one  cycle  once  printing  began. 
The  magnitude  of  these  peaks  were  not  always  identical.   This  indicated 
that  the  system  had  not  reached  steady  state;  however,  in  most  cases, 
the  difference  between  the  peaks  was  less  than  plotting  accuracy.  The 


average  of  the  magnitude  of  these  peaks  was  taken  to  be 


a 


With  values  of  the  damping  ratio  greater  than  about  .4,  the  method 


of  obtaining 


a 


as  described  above  was  used.   At  damping  ratios  of 
about  0.3  and  lower,  there  began  to  appear  an  oscillation  at  the  natural 
frequency  of  the  system  super- imposed  on  the  output  wave  form.  Its 
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presence  could  be  expected  when  a  limit  cycle  would  be  predicted  by 
Ref.  (2). 

Fig.  k- la  shows  this  phenomenon.   The  damping  ratio,  u     ,  was 
0.1  in  this  case.   Shown  in  Fig.  ^4-- lb  is  the  output  waveform  of  the 
same  system  with  a  higher  damping  ratio.   (  ^f    =  0.5)   It  can  be 
seen  here  that  there  is  no  super- imposed   limit  cycle  .   To  obtain 

0C   for  the  system  whose  outpuf.  is  shown  in  Fig.  ^- la,  it  was  nec- 
essary to  take  the  mean  of  the  super- imposed  oscillation.   In  the  case 
shown,  this  was  not  particularly  difficult.  However,  as  co  was  in- 
creased towards  the  undamped  natural  frequency,  this  procedure  became 
more  difficult.   Fig.  i+-2a  shows  the  output  waveform  for  the  same 
system  at  a)  =  .6.   From  this  figure,  it  can  be  seen  that  it  is  vir- 
tually impossible  to  establish  the  mean  of  the   limit  cycle   oscilla- 
tions.  Even  if  the  output  wave  shape  were  available  for  a  longer  time 
period,  it  would  still  be  virtually  impossible  to  obtain  an  accurate 
mean  value.   To  compare  the  output  at  the  same  forcing  frequency,  but 
with  higher  damping  ratio,  Fig,  *4—  2b  was  drawn.  Here  y   =0.5>  and  it 
can  be  seen  that  there  is  no  difficulty  picking  the  peak  output  value. 

Since  the  method  of  picking  |  QA     varied,  it  was  not  possible  to 


program  the  selection  of   C7C 


In  some  cases,  a  programmed  method 


would  have  been  satisfactory,  but  unless  all  the  data  were  examined, 


the  value  of 


fc)c   thus  obtained  would  be  questionable. 
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k  .2  Frequency  Response  Variation  with  Damping  Ratio 

Frequency  response  curves,  for  plastic  impact  of  the  gear  teeth, 

were  obtained  for  values  of-^-of  10,  5  and  0.95>  with  damping  ratio 

Jf  varying  from  0.1  to  0,8.   These  curves  are  shown  in  Figures  4- 3> 

h-k   and  4-5  respectively.   For  comparison  purposes,  a  similar  family 

of  curves  for  a  comparable  linear  system  is  shown  in  Fig.  k-6. 

It  was  determined  by  test  runs,  where  both  A  and  A  were  varied 

proportionally,  that  the  frequency  response  was  identical  as  long  as 

the  ratio  —  was  cdnstant.   For  example,  the  same  frequency  response 

resulted  with  A  =  1.0  and  A  =  0.1  as  with  A  =  2.0  and  A  m   0.2.   This 

makes  the  frequency  response  results  generally  applicable  to  all 

simple  second  order  servo  systems  for  the  same  relative  amounts  of 

backlash  since  the  system  was  normalized  with  cc  =1.   For  this  reason 

n 

the  abscissae  of  the  curves  are  labeled  ayu)  ,  and  the  curves  can  be 

n' 

applied  to  any  simple  second  order  system  by  multiplying  by  a)  of  the 
system . 

The  resonant  frequencies,  for  values  of  damping  ratio  that  resulted 
in  resonance  ( J  ^  0.6),  were  found  to  increase  as  the  amount  of  back- 
lash was  increased.   This  increase  in  resonant  frequency  was  greatest 
for  the  higher  values  of  ^   while  it  was  almost  insignificant  at  a 

y  of  o.i. 

The  peak  values  of  the  ratio  of  output  position  to  input  ampli- 
tude at  resonance  M  ,  increased  above  the  linear  peaks  as  the  amount 

P 

of  backlash  was  increased.   The  increase  was  greatest  at  the  lowest 
value  of  y  .   A  plot  showing  this  variation  from  the  linear  for  dif- 
ferent amounts  of  backlash  is  shown  in  Fig.  4- 7» 
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As  might  be  expected,  the  effect  of  backlash  is  greatest  at  the 
resonance  frequencies.   The  frequency  response  curves,  for  all  amounts 
of  backlash,  are  above  the  linear  curve  at  values  of  a>  below  the  re- 
sonant frequencies.  At  frequencies  greater  than  the  undamped  natural 
frequency  and  for  values  of  !/)  0.4,  the  response  with  backlash  drops 
off  more  rapidly  than  the  linear  and  the  degree  of  drop-off  is  greater 
as  the  backlash  is  increased.   In  the  same  frequency  range,  the  response 
curve  at  low  values  of  tf    is  slightly  above  the  linear  curve  (ignoring 
sub-harmonic  resonance  which  will  be  discussed  later).   This  indicates 
that  with  heavy  damping,  backlash  results  in  a  more  rapid  attenuation 
of  output  at  frequencies  above  the  undamped  natural  frequency  of  the 
system,  while  with  light  damping,  attenuation  is  slightly  less  than  Lhe 

linear  up  to  2  (i)  . 
n 
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h . 3   E ffect  of  Magnitude  of  Backlash  on  Frequency  Response  for  Constant  J 

A  more  detailed  examination  of  the  effect  of  the  amount  of  backlash 
in  the  system  compared  with  the  input  amplitude  was  made  for  a  constant 
damping  ratio.   The  value  of  ^f   =  0.5  was  chosen  since  it  represents  a 
practical  damping  ratio  which  is  less  than  optimum  and  results  in  a 
moderately  oscillatory  system.  Also,  no  limit  cycle  exists  at  this 
value  of  y   according  to  Anderson  and  Luckett,  Ref.  2.   It  was  pointed 
out  earlier  that  as  long  as  the  ratio  of  A/^  is  constant  the  frequency 
response  is  invariant.   Therefore,  A/\   is  a  logical  system  parameter. 

Figure  U-8  shows  the  frequency  responses  at  ^f   =  0.5  for  r\/£±     of 
10,  5  and  0.95  with  plastic  impact  of  the  gear  teeth.   The  linear  re- 
sponse is  also  shown  on  the  plot.   It  can  be  seen  that  the  response 
deviated  from  the  linear  even  when  the  magnitude  of  backlash  was  only 
one- tenth  of  the  input  amplitude.   The  deviation  increased  as  the  rela- 
tive magnitude  of  backlash  was  increased.   The  effect  of  backlash  was 
to  change  the  response  curve  to  resemble  the  response  for  a  lower  damp- 
ing ratio. 

At  frequencies  above  the  undamped  natural  frequency,  with 
A  /A  =  0.95>  the  response  curve  dropped  below  the  linear  curve  and 
reached  a  minimum  near  a>  =  2 .   Up  to  this  input  frequency  the  ratio 
of  output  frequency  to  input  frequency  was  one  to  one.  With  a  slightly 
higher  cd,  the  ratio  of  output  to  input  frequencies  f  /f .  rapidly  tran- 
sitioned to  1/2  and  the  magnitude  of  the  output  abruptly  increased. 
This  was  the  first  of  a  series  of  sub-harmonic  resonance  peaks  which 
were  encountered  as  co  was  increased  further.   The  actual  values  of  M 
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between  sub-harmonic  resonance  peaks  were  not  all  definitive  as  will  be 
shown  in  the  discussion  of  the  output  wave  shapes  for  this  region  of  a). 
At  some  frequencies  where  input  frequency  was  not  an  integer  multiple 
of  the  output  frequency,  the  output  magnitude  could  not  be  defined 
because  of  the  random,  unsteady  nature  of  the  load  position.   This  is 
shown  by  a  dashed  line  in  Fig.  4-8. 

4 .4   Sub-harmonic  Resonance  with  Plastic  Impact 

The  first  sub-harmonic  resonance  peak  occurred  at  go  =  2.2,  (Fig.  k-Q). 
At  input  frequencies  of  2.0  and  below,  the  outputs  were  steady  with  ampli- 
tudes decreasing  with  increasing  co,  and  the  ratios  of  output  to  input 
frequencies  were  l/l .   The  wave  shapes  of  the  output  for  oj  =  1.8>  1.9 
and  2.0  are  shown  in  Figs.  4-9  and  4-10.  At  a>  =  2 . 1  the  output  was 
very  unsteady  (Fig.  4- 10b )  and  a  value  for  M  could  not  be  determined. 
At  this  input  frequency,  the  ratio  of  output  frequency  to  input  frequency 

f  It.   was  not  constant  but  fluctuated  with  an  average  of  between  1/1  and 
o   i 

1/2. 

When  oo  was  increased  to  2.2,  the  frequency  ratio  f  /f .  became  a 
steady  l/2.   The  output  wave  shape  became  steady  but  unsymmetrical,  as 
shown  in  Fig.  4- 11a.   The  peaks  of  this  output  wave  occurred  at  the  times 
of  contact  of  the  gear  teeth.   The  reason  for  the  unsymmetrical  wave 
shape  was  due  to  the  difference  in  relative  velocities  of  the  motor  and 
load  at  the  time  of  contact.   Referring  to  Fig.  4- 11a,  the  higher 
(positive)  peaks  occurred  as  a  result  of  the  lower  relative  velocity 
at  the  instant  prior  to  the  preceding  contact  (smaller  peaks).   This 
smaller  relative  velocity  imparted  a  higher  velocity  to  the  load  which 
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caused  a  greater  excursion.  The  opposite  was  true  in  going  from  contact 

at  the  higher  peaks  (high  relative  velocity)  to  the  smaller  peaks .  This 

phenomenon  is  believed  to  be  related  to  the  fact  that  there  was  plastic 

impact  of  the  gear  teeth.   The  value  of  M  at  the  sub- harmonic  resonant 

peak  at  {jd  =  2.2  in  Fig.  4-8  came  from  the  average  of  the  high  and  low 

peaks . 

When  (D  was  increased  to  2.3  the  output  wave  became  erratic,  the 

ratio  f  It.    was  transitioning  from  l/2  to  l/3  and  the  exact  value  of  M 
01 

was  indeterminant  but  was  obviously  lest  than  at  CO  =2.2  (Fig.  4- lib). 

At  go  =  2.4  the  output  wave  (Fig.  4~l2a)  was  similar  to  the  wave 
shape  at  the  previous  sub-harmonic  resonance  peak  of  co  =  2.2.   The  wave 
was  steady  but  unsymmetrical  and  the  ratio  of  f  /f,  was  a  fixed  1/3. 
The  value  of  M  was  further  reduced  (averaged  from  peaks  of  output  wave). 

Going  to  a)  =  2.5  (Fig.  4- 12b)  the  wave  shape  was  again  unsteady 

but  f  /f.  was  a  consistent  1/3  and  M  began  to  increase.   This  increase 
01 

in  M  continued  for  the  next  four  a)  s. 

Wave  shapes  for  cu' s  2.6,  2.7>  2.8  and  2.9  shown  in  Figs.  4-13  and 
4-14,  show  very  convincingly  a  sub-harmonic  resonance.  At  u)  =  2,6,  the 
general  wave  shape  of  the  output  during  this  resonance  is  first  observed 
For  all  four  of  these  frequencies,  the  system  output  had  the  same  wave 
shape.  As  the  frequency  was  increased,  the  M  value  also  increased. 
All  four  of  these  frequencies  gave  a  f  /f.  of  l/3«   Just  as  the  output 
wave  shape  was  very  steady,  the  f  /f,  was  equally  as  steady  at  l/3° 
The  sub-harmonic  rise  associated  with  these  frequencies  is  shown  on 
Fig.  4-8.   Since  the  value  of  M  for  these  frequencies  was  so  well 
defined,  this  portion  of  the  curve  was  not  dashed  as  was  required  for 
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frequencies  between  a)  =  2.0  and  co  -  2  .4  . 

A  further  increase  in  frequency  beyond  co  -  2.9  resulted  in  a  drop- 
off in  M.   Fig.  4-15  shows  the  output  wave  shape  for  co  =  3-0.   This 
wave  shape  was  taken  for  a  period  of  over  50  time  constants.  As  can  be 
seen,  there  are  portions  of  this  wave  which  are  quite  steady  and 
portions  which  are  very  unsteady.   The  steady  portion  had  a  f  /f.  of 
l/i.   The  M  value  taken  from  this  steady  portion  was  used  to  plot  on 
Fig.  4-8.   No  M  value  could  be  obtained  from  the  unsteady  portion  due 
to  its  inconsistent  movement.   The  f  /f.  during  the  unsteady  portions 
was  not  defined  either.   It  averaged  between  1/3  to  l/4 , 

At  frequencies  greater  than  3-0,  the  output  motion  was  erratic. 
This  can  be  seen  in  Fig.  4-16  for  co  =  3»1  and  4,0  as  well  as  in 
Fig.  4-17  for  co  =  5«0.   Since  M  could  not  be  deduced  from  these  wave 
shapes,  no  points  are  plotted  on  Fig.  4-8  beyond  co  =  3°0c 

The  f  /f .  for  co  =  3-l>  Fig.  4-l6a,  was  l/3«   This  may  be  mislead- 
ing because  only  a  small  amount  of  the  output  wave  was  obtained.   At 
CO  =  4.0,  Fig.  4- 16b,  the  f  /f.  was  between  l/k   and  1/5  .   Fig.  4-17, 
for  co  =  5^0,  shows  the  output  wave  shape  for  a  considerable  time.   It 

can  be  seen  here  that  the  wave  was  not  consistent.   The  f  /f.  for  this 

o      x 

co  varied  widely.   Depending  on  the  portion  of  the  curve  examined,  the 
f  /f.  was  as  small  as  l/6  and  as  large  as  l/4. 

Although  no  sub-harmonic  resonance  was  found  beyond  co  =  2 /„->,  it 
would  appear  that  it  does  exist.   In  transitioning  from  one  sub- 
harmonic  frequency  to  another,  the  output  wave  was  quite  unsteady. 
When  a  definite  output  frequency  was  established,  the  M  was  also 
defined.   Since  the  f  /f .  ratio  changed  from  1/3  at  co  =  2-9  to  l/6  at 
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CD  =  5 «0>  a  sub-harmonic  resonance  at  l/4  and  l/5  must  have  occurred  between 
these  two  frequencies. 

From  the  limited  data  available,  it  appears  that  the  system  tends 
to  resonate  at  odd  sub-harmonics .   Frequencies  between  co  =  2.0  and 

O)  =  2.2  did  not  have  a  fixed  f  /f .  .   Then  at  co  =  2.2,  a  steady  l/2  ratio 

01  ' 

was  observed. 

During  the  build  up  to  the  third  sub-harmonic  resonance  peak,  a 
steady  ratio  of  l/3  was  present.   This  steady  f  /f.  over  several  values 
of  frequency,  leads  one  to  believe  the  system  would  tend  to  seek  an  odd 
sub-harmonic  frequency.   Immediately  past  the  sub-harmonic  peaks,  the 
system  began  seeking  the  next  sub-harmonic  frequency.   The  even  sub- 
harmonics  did  not  appear  steady  except  at  co  =  2,2.   There  was  a  con- 
siderable range,  co  -   2,k   to  a)  =  2.9>  where  f  ft.   was  a  steady  1/3  add- 
ing support  to  the  belief  that  the  system  would  seek  odd  sub-harmonics. 
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k  .5   Effect  of  Elastic  Impac 

Frequency  response  data  were,  obtained,  at  a  constant  f  =  0.^,  for 
values  of  A/A  equal  to  10,  5  and  0.95  with  e  =  0.6  and  e  -  0.8.   These 
parameters  were  the  same  as  for  the  results  discussed  in  Section  h.3f 
except  for  the  values  of  the  coefficient  of  restitution  e.   The  response 
curves  for  e  -  0.6  and  e  =  0.8  are  shown  in  Figs.  k~ 18  and  U-I9.   The 
values  of  M  for  all  three  values  of  A/A  were  found  to  be  slightly 
less  for  both  cases  of  elastic  impact  than  for  e  =  0  (Fig.  h~Q)   at 
frequencies  near  resonance.   This  indicates  that  elastic  impact  tends 
to  reduce  M  as  compared  with  plastic  impact.   There  was  practically  no 
difference  between  the  response  curves  for  e  =  0.6  and  e  -  0.8  up  to 
input  frequencies  of  03  =  2. 

With  a/A  =  0.95  sub-harmonic  resonance  began  to  occur  at  about 
CO  =  2.8  for  e  =  0.6.   Up  to  0)  =  2.7  (Fig,  4-20a)  the  ratio  f  /f.  was 
a  steady  1/1  and  the  output  wave  shapes  were  steady  and  symmetrical. 
At  a)  =  2.8  the  output  waveshape  became  erratic  (Fig.  U-21a)  and  f  /f . 
was  transitioning  from  1/1  to  l/2.   The  response  curve  remained  above 
the  linear  at  all  frequencies  (Fig.  4-18),  unlike  the  equivalent  curve 
for  e  =  0  (Fig.  h-Q)   which  dropped  below  the  linear  between  M  and  the 
occurrence  of  sub-harmonic  resonance. 

The  response  curve  reached  a  minimum  at  O)  =  2.8.  At  the  next 
test  frequency  a)  =  2.9>  the  output  increased  and  the  ratio  f  /f.  was 
between  l/2  and  1/3.   No  even  harmonic  was  found  as  there  was  for  e  =  0. 
The  wave  shape  at  O)  =  2.9  was  also  uneven  with  very  inconsistent  peaks 
(Fig.  4-20b).  At  0)  =  3.0  and  a)  =  3.2,  f  /f.  was  still  between  l/2  and 
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1/3  the  value  of  M  was  increasing  with  unsteady  waveshapes  (Figs,  k- gib 

and  4-22a).  At  a)  -   J>M   (Fig.  4-22b)  f  /f .  became  a  consistent  l/3  and 

M  remained  about  the  same,  but  the  wave  shape  was  still  not  steady.  As 

O)  was  increased  to  3»7  an<3  to  k  .0   (Fig,  4-23)  the  waves  became  symmetrical, 

f  /f .  remained  l/3  and  the  value  of  M  was  decreasing.  At  the  highest 
o   1 

test  frequency  a)  =  4.3,  M  was  still  decreasing,  the  wave  shape  was  be- 
coming unsymmetrical  (Fig.  4-24)  and  f  /f.  was  decreasing  toward  l/4 . 

•    o   i 

In  comparing  the  response  curves  of  Fig.  4-8  (e  =  0  and  Fig.  4-18 
(e  =  0.6)  for  a/a  =  0.95*  it  can  be  seen  that  the  addition  of  bounce 
changed  the  curve  and  the  sub-harmonic  resonance  pattern  at  frequencies 
above  0)  .   Bounce  eliminated  the  drop  of  the  response  curve  below  the 
linear  response  and  apparently  al  o  eliminated  the   even  sub-harmonic 
resonance  peak.   The  odd  (l/3)  sub-harmonic  resonance  peak  was  lower 
with  bounce  but  was  wider,  i.e.  the  value  of  f  /f.  =  l/3  persisted  over 
a  wider  range  of  frequencies.   Sub-harmonic  resonance  started  at  a 
higher  frequency  when  bounce  was  introduced. 

For  e  =  0.8  and  a/a  =  0.95>  sub-harmonic  resonance  started  to 
occur  at  about  the  same  value  of  a)  as  for  e  =  0.6.   The  investigation 
did  not  go  beyond  this  frequency  for  e  =  0.8. 

In  general,  elastic  impact  of  the  gear  teeth  tended  to  reduce  the 
effects  of  the  nonlinearity  for  both  the  magnitude  of  the  output  and  the 
amount  of  sub-harmonic  resonance. 

4 .6  Response  Characteristics  with  Very  Low  Damping  Ratio 

In  an  effort  to  find  jump  resonance  and  the  factors  influencing 
its  occurrence,  some  runs  were  made  with  a  damping  ratio,  y   ,  equal  to 
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0,05.   Considerable  computer  time  was  required  for  calculating  the  out- 
put at  such  a  low  damping  ratio.   Therefore,  only  a  limited  number  of 
frequencies  were  investigated. 

The  basic  parameters  for  this  investigation  were:   damping  ratio 
equal  to  0.05>  a/a  equal  to  1.2,  coefficient  of  restitution  equal  to 
zero,   The  frequencies  used  were  between  0.5  and  1.8  with  emphasis  on 
the  frequencies  near  the  undamped  natural  frequency.   Results  of  these 
runs  did  not  show  any  jump  resonance;  however,  some  interesting  things 
were  shown. 

For  a  linear  system,  the  resonant  frequency  is  always  less  than  the 
undamped  natural  frequency.   (Resonant  frequency  is  the  frequency  at 
which  the  resonant  peak  occurs  on  a  frequency  response.)  Regardless 
of  how  small  the  damping  ratio  gets,  the  resonant  frequency  is  always 
less  than  the  undamped  natural  frequency;  however,  as  y    approaches 
zero  the  resonant  frequency  approaches  but  never  exceeds  the  undamped 
natural  frequency. 

With  backlash  present,  this  is  not  always  so.   Fig.  4-25  shows  a 
resonant  peak  at  a  frequency  between  1 .01  and  1.02  which  was  greater 
than  0)  o   it  has  been  shown  previously  that  the  effect  of  increasing 
backlash  was  similar  to  reducing  the  damping  ratio.   In  addition  to 
this  effect,  it  appears  that  with  a  sufficiently  large  backlash  and 
small  enough  jf  ,  the  resonant  frequency  is  increased  more  than  the  in- 
crease associated  with  a  reduction  in  ^f    alone.   This  may  not  be 
limited  to  low  damping  ratios.   It  may  have  occurred  at  higher  damping 
ratios  but  was  masked  due  to  insufficient  data  being  taken  in  the 
vicinity  of  the  resonant  peak. 
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The  height  of  the  resonant  peak  for  this  low  ^f  is  considerably 
higher  than  for  the  linear  case.   However,  extrapolation  of  Fig.  4-7 
gives  approximately  the  same  value  for  this  peak  as  shown  on  Fig.  4-25. 
At  frequencies  near  the  resonant  frequency,  between  0,95  a"d  1.05,  the 
output  wave  shapes  were  steady .   The  frequency  ratio  was  also  steady  at 
one  to  one.   Fig.  4-26  shows  the  output  wave  shape  at  co  =  0.95  and 
Fig.  4-2?  shows  it  at  CD   -  1.0.  As  *can  be  seen,  these  are  steady. 

The  output  below  a>  -  0„9  was  rather  difficult  to  interpret.   By 
the  previously  discussed  methods  for  obtaining  M,  it  was  not  possible 
to  determine  this  quantity  accurately.   Therefore,  values  of  M  for  these 
frequencies  were  not  plotted  on  Fig.  4-25. 

Fig.  4-28  shows  the  wave  shape  for  O)  -  1.2.   Even  though  the  output 
wave  shape  is  steady  at  this  frequency,  the  input  frequency  is  not  an 
even  multiple  of  the  output  frequency.   This  uneven  relationship  has 
been  mentioned  previously  (Section  4.4)  as  an  indication  of  the  onset 
of  sub-harmonic  resonance.   Therefore,  the  low  damping  ratio  decreases 
the  frequency  at  which  sub-harmonic  resonance  appears.   Insufficient  data 
are  available  to  determine  the  exact  location  of  the  first  sub-harmonic 
resonant  peak,  but  it  appears  to  occur  at  a  frequency  near  u)  =  1.6.   The 
M  at  this  peak  is  approximately  9?  and  the  frequenc)  ratio  is  even  at 
one  to  two . 
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If  .7  Test  of  Program  and  Validity  of  Results 

For  there  to  be  complete  confidence  in  the  results  of  this  investi- 
gation, some  mention  should  be  made  of  the  steps  taken  to  prove  the  com- 
puter program  from  which  the  results  were  obtained.  Also,  some  mention 
of  the  validity  of  the  results  is  in  order.   It  is  the  purpose  of  this 
section  to  discuss  these  subjects. 

Since  there  was  no  easy  way  that  the  system  investigated  could  be 
simulated  on  an  analog  computer;  and  since  little  information  is  avail- 
able on  such  non- linear  systems,  it  was  necessary  to  test  the  program 
by  approximation.   First  of  all,  the  program  was  divided  into  two  basic 
parts.   These  were,  the  combined  system  portion  and  the  separated  system 
portion.  As  integral  parts  of  these  two  sections  were  the  basic  decisions 
which  had  to  be  made  to  transfer  from  one  to  the  other. 

It  was  rather  simple  to  test  the  combined  section.   This  was  done 
by  setting  the  backlash,  A  ,  equal  to  zero.   Doing  this  effectively 
eliminated  the  separated  system  portion  and  the  decisions  required  to 
transfer  into  that  portion.  With  the  backlash  equal  to  zero,  the  output 
should  be  the  same  as  the  output  from  a  linear  system.   Comparison  of 
the  M  obtained  from  the  print  out  data,  as  described  in  section  ^—1, 
with  linear  results  obtained  by  conventional  means  were  very  favorable. 
The  slight  differences  could  be  attributed  to  the  computer  solution 
still  being  in  a  transient  state.   The  data  from  the  computer  solution 
was  taken  after  five  time  constants  had  elapsed.   The  following  table 
shows  this  comparison. 
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TABLE  IV- 1 


tf«  .8 


A,  o 


03 

M,  computer 

M,  linear 

.1 

.997^ 

.99716 

.2 

.98810 

.98821 

•  3 

.97172 

.97197 

,k 

=9^693 

.9^69^ 

•  5 

,91180 

.91192 

It  was  not  possible  to  test  the  separated  portion  alone.   There- 
fore, it  was  necessary  to  test  this  portion  as  part  of  the  overall  test 
of  the  program.   For  a  linear  system,  the  ratio  of  forcing  function  am- 
plitude to  the  amount  of  backlash,  A/A  ,  is  infinite  since  the  backlash 
is  zero.   Therefore,  if  this  infinite  ratio  is  approximated  by  a  non- 
linear system,  the  results  should  be  very  close  to  linear  results.   To 
approximate  the  linear  system,  an  A,/A  of  100  was  chosen.   The  results 
of  this  test  run  were  within  2.5  percent  of  the  linear  response  at  all 
frequencies  less  than  the  resonance  frequency.  At  frequencies  above 
resonance,  the  non- linear  response  was  less  than  five  percent  different 
than  the  linear  response  out  to  the  highest  frequency  used.   The  maximum 
percent  difference  occurred  at  the  highest  frequencies  where  small  dif- 
ferences account  for  higher  percentages. 

Since  the  program  was  capable  of  approximating  a  linear  system 
while  actually  going  through  all  of  the  calculations  necessary  for  a  non- 
linear system;  it  was  reasonable  to  assume  that  with  a  greater  amount 
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of  backlash,  smaller  — — ,  the  program  would  give  valid  results. 

Even  though  the  time  solution  of  the  differential  equations  was 
considered  accurate,  the  method  of  determining  M  may  have  introduced 
inaccuracies.  For  damping  ratios  greater  than  0.3,  the  method  of 
obtaining  M  was  quite  accurate.   The  only  improvement  would  have  been 
to  calculate  for  a  longer  time,  but  this  would  have  required  more  time 
than  could  be  justified  by  the  increased  accuracy. 

At  damping  ratios  less  than  0.3?  the  interpretation  of  the  output 

waveform,  kjc   ,  was  not  so  clear  cut.  At  these  damping  ratios,  the 

ii,  .  ,      ,  it  .      ,     , 

limit  cycle  was  super- imposed  on  the  output  wave.   At  some  frequencies, 

O)  less  than  about  0A?  the  method  used  to  obtain  M  was  quite  good. 

(Mean  of   limit  cycle  oscillations  )  However,  for  frequencies  between 

about  0.^4-  and  0.7>  the  M  was  not  so  nicely  defined.   It  was  in  this 

frequency  range  for  low  if  S  that  the  results  were  not  as  accurate  as 

they  were  otherwise.  At  all  frequencies  above  an  a)  of  about  0.7>  the 

results  wert  quite  good  except  as  dicussed  in  the  sub-harmonic  section. 

Mention  has  been  made  that  data  was  taken  after  five  time  constants 

based  on  a  linear  system.   The  wave  shapes  at  earlier  times  were  looked 

at  to  determine  if  this  delay  were  required*   From  these,  it  was  evident 

that  the  system  was  in  a  transient  state .   Even  after  five  time  constants, 

there  was  evidence  of  transients.  However,  the  transient  affect  was 

quite  small.   Depending  upon  the  parameters,  the  transient  was  virtually 

gone  shortly  after  five  time  constants  had  elapsed.  The  only  exception 

to  this  was  discussed  as  part  of  sub-harmonic  resonance.   In  the  cases 

where  steady  state  was  reached,  it  was  reached  in  less  than  six  time 

constants . 
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k  ,8  Application  of  Results  and  Extension  of  Method 

While  there  is  limited  direct  application  of  the  results  of  this 
investigation  to  servo  design,  there  are  some  points  that  could  prove 
helpful  in  the  design  of  certain  types  of  systems.   If  a  system  with 
backlash  were  to  be  operated  with  a  periodic  input,  Fig.  U~7  would  be 
useful  in  determining  how  much  backlash  could  be  tolerated  for  a  given 
input  amplitude  to  meet  specified  peak  outputs.   Also,  the  results 
showing  the  effect  of  gear  resilience  can  be  of  benefit  in  choosing 
a  gear  material,  since  higher  resilience  tends  to  reduce  the  effect  of 
the  backlash. 

The  general  form  of  the  computer  program  can  be  adapted  to  calcu- 
late any  particular  frequency  response  of  similar  second  order  systems 
with  known  system  constants  (K,  J  ,  J  ,  A,  etc.)  for  specified  para- 
meters.  The  program  can  be  very  easily  modified  to  obtain  transient 
response  of  second  order  systems  with  any  desired  sets  of  constants  and 
parameters.   Similarly,  it  can  be  readily  applied  to  steady  state  (limit 
cycle)  investigations. 

As  has  been  pointed  out,  the  method  of  obtaining  frequency  response 
results  in  this  investigation  requires  a  great  deal  of  digital  computer 
time.   If  further  investigations  were  to  be  based  on  the  same  general 
approach,  especially  at  very  low  damping  ratios,  a  technique  of  going 
from  one  frequency  to  another  without  having  to  go  through  a  new  tran- 
sient period  each  time,  should  be  devised.   Such  a  program  was  briefly 
attempted,  as  mentioned  in  Section  3>  which  transferred  boundary  con- 
ditions of  the  differential  equations  at  the  first  frequency  to  the 


appropriate  differential  equations  as  initial  conditions  at  the  next 

ii      ii 
frequency.   This  was  an  attempt  to  utilize  the   energy  stored  in  the 

system  to  more  rapidly  set  up  a  steady  state  condition  for  the  next 

frequency.   It  appears  this  concept  may  be  the  only  way  that  would  make 

feasible  this  digital  method  to  more  extensive  investigations.   It  could 

result  in  showing  jump  resonance  when  it  exists,  where  the  program  used 

herein  apparently  did  not. 
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Fig.   4-13 
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Fig.   4-23 
OUTPUT  WAVE  SHAPES 
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5.  CONCLUSIONS 

Based  on  the  findings  of  this  investigation,  the  following  conclu- 
sions are  made: 

1.  The  system  frequency  response  remains  unchanged  if  A  and  A 
are  varied  proportionally,  with  all  other  parameters  constant . 

2.  The  presence  of  backlash  in  the  system  has  an  effect  similar 
to  lowering  the  damping  ratio  of  the  system. 

3.  The  effect  of  backlash  is  more  pronounced  at  lower  damping 
ratios. 

k.      Increasing  the  backlash  causes  a  more  oscillatory  response. 

5  „  With  a  very  low  damping  ratio  and  a  small  value  of  A/A  ,  the 
first  resonant  peak  can  occur  at  a  frequency  greater  than  the  undamped 
natural  frequency. 

6.  At  sub-harmonic  resonance  the  ratio  of  input  frequency  to  output 
frequency  is  an  integer. 

7.  The  output  magnitude  is  erratic  when  the  input  frequency  is 
not  an  integer  multiple  of  the  output  frequency c 

8.  The  odd  sub-harmonics  are  more  pronounced;  and  with  elastic 
impact,  the  even  sub-harmonic  does  not  appear, 

9*  The  lower  the  damping  ratio,  the  lower  the  frequency  at  which 
sub-harmonic  resonance  will  occur. 

10.  With  all  other  parameters  constant,  sub-harmonic  resonance 
will  occur  at  a  lower  frequency  with  plastic  impact  than  with  elastic 
impact. 
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11.  The  effect  of  elastic  impact  is  to  reduce  the  effect  of  back- 
lash. 

12.  With  a  coefficient  of  restitution  of  either  0.6  or  0.8,  the 
frequency  response  of  the  system  is  essentially  the  same. 

13.  Combinations  of  low  damping  ratio  and  high  frequency  require 
excessive  computer  time  to  obtain  the  frequency  response. 

lU .  When  a  limit  cycle  exists  (Ref.  2),  an  oscillation  at  the 
undamped  natural  frequency  of  the  system  is  superimposed  on  the  output 
motion. 
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APPENDIX  A  -  FLOW  CHART  OF  COMPUTER  OPERATION 
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EXPLANATION  OF  FLOW  CHART  NUMBERS 

I.   Set  initial  conditions  and  assign  values  to  constants  and  parameters 
o£  system  equations:   i,e„,  J,     |M  ,   fL,  A,  A  ,  es  K.,  jO   ,  cd  , 
Problem  time  for  delay  of  live  time  constaatSj,  number  of  input 
cycles  desired  after  five  time  constants, 

2*  *Print  heading  for  results, 

3.  *Print  run  number  and  values  of  e,  jf  and  A. 

4.  *Print  column  headings  for  Time,  0T ,  0T  ,  0%,„  0^  and  C?  . 

L7    L    M    M       R 

5.  Set  first  initial  conditions  for  combined  differential  equation- 

* 

Time  =  0,  ©  =  0,  Q   =0 

6.  Select  values  of  o>  and  At 

7»  *Print  values  of  u),  At,  and  J  . 

8.   Solve  differential  equation  for  combined  operation 

9»   Increment  time  by  integration  interval   t  — »-  t  +  At  „ 

10.  Calculate  value  of  input  0  . 

R 

11.  Does  time  equal  five  time  constants'? 

12.  Is  backlash  zero? 

13.  Solve  for 


ec  =  -2yconec  -^n6c  +  ujn  A  sin  cot 


14.      Calculate 


N  = 


-  QcA 


6 


C 


15.  Calculate  /     / 

16.  Does    Nlc  =   -  NL? 
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17.  Count  a  time  increment  At, 

18.  Gone  through  IOC  At6 s? 

19.  Set  At  count  to  zero. 

20.  *Print  values  of  Time,  0  .  9  and  0„ „ 

*    c'        c       R 

21.  Does  time  equal  five  time  constants  plus  time  for  desired  number 
of  input  cycles  after  transient  period1? 

22.  Set  initial  conditions  for  second  and  third  differential  equations. 

23.  Does  time  equal  five  time  constants? 

2k.     *Print  values  of  Time,  0  ,  0  and  0„ . 

'    c'    c       R 

25.  Are  gear  teeth  separating  from  original  contact  position  plus  A» 

26.  Calculate     0   =   ®L " ^ 

27.  Calculate     A  ~     Ji- 

28.  Solve  differential  equation  for  motor  alone 

vJm       Um       v^m 

29.  Solve  differential  equation  for  load  alone 

30.  Increment  time   t  -*•  t  +  At, 

31.  Calculate  value  of  input  ©  , 

32.  Does  time  equal  five  time  constants? 

33.  Count  a  time  increment  At. 
jU.  Gone  through  100  At'sV 
35«  Set  At  count  to  zero. 

36.  *Print  values  of  Time,  0  ,  0  ,  0  ,  0  and  Q    . 
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37.  Does  time  equal  five  time  constants  plus  time  for  desired  number 
of  input  cycles  after  transient  period? 

38.  Are  gear  teeth  touching  in  initial  contact  position? 

39.  Are  gear  teeth  touching  in  initial  position  plus  A? 
kO.      Calculate  load  bounce  velocity  Q. 

Jm 


•  / 


eL  = 


JmV'Jl 


^eM(,+  e)  +  9L^-e) 


J 


Ul»   Calculate  motor  bounce  velocity  Q 

A'  -    6l  -e/>9M  +  e  6L 


k2 ,  Do  gear  teeth  stay  together? 

^4-3.  Set  initial  conditions  for  combined  differential  equation 

kk.  Does  time  equal  five  time  constants? 

45.  *Print  values  of  Time^  8  ,  ©L#  0,  B  and  ©  „ 

k6.  Go  to  combined  differential  equation. 

kj .  Completed  calculation  for  all  values  of  a)? 

W3.  Stop. 


i>     11 
*   Print   refers  to  putting  the  indicated  items  or  their  values  onto 

magnetic  tape  for  later  printing. 
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APPENDIX  B 
COMPUTER  FORTRAN  SOURCE  PROGRAM 


PROGRAM  BLASH 

DIMENSION  DM(3),  DL(3),  W(23),  TINT128).  Y1I2),  Y2 ( 2  ) ,  Y3(2) 
COMMON  ZETA,  WN,  A,  OMEGA,  FLf  DLOAD,  FM,  DMOTORf  GC,  THEL 
ZETA  =  0.5 


tfc 


:1 


DELAY  =  10.0 
A  =  1.0 
BL  =  0.  I 
E  =  0.0 
CYC  =  6.3 
GC  =  1.0 
RHO  =  1.0 
WN  =  1.0 

EPS  =  .45474735^-12 
1   FORMAT  (3F10.I) 

READ  INPUT  TAPE  7,  I  ( DL ( J )  ,  J  =  1,  3) 
READ  INPUT  TAPE  7,  I  { DM ( J ) ,  J  =  1,  3) 

23  FORMAT  (7F10.2) 

READ  INPUT  TAPE  7,  23  (W(I),  I  =  1,  28) 

24  FORMAT  (6F12.10) 

READ  INPUT  TAPE  7,  24  (MNT(I)t  I  =  1,  28) 
50    FORMAT  (58H1RESULTS  OF  NONLINEAR  SERVO  WITH  BACKLASH    GRAHAM  /  LL 
10YD) 
WRITE  OUTPUT  TAPE  8,  50 
2   OFORMAT  ( 23H0PR00UCTI0N  RUN  19  E  =  , 

IF3.1,  3X,  7HZETA  =  ,  F4.2,  3X,  4HA  =  ,  F6.3) 
WRITE  OUTPUT  TAPE  8,  2,  E,  ZETA,  A 

3  OFORMAT  (5X,  UHTIME,  12X,  9H0UT  POSIT,  13X,  7H0UT  VEL,  9X,  1IHM0T0R 
1  POSIT,  1JX,  9HM0T0R  VEL,  11X,  9HINPUT  POS) 

WRITE  OUTPUT  TAPE  8,  3 
DO  19  J  =  2,2 
DO  19  I  =  8,  9 
PAUSE 

DELT  =  TINT( I) 
OMEGA  =  W( I) 
DLOAD  =  DL(J) 
DMOTOR  =  DM(J) 
FINI  =  CYC/WII ) 
T  =  0.0 
Yl(  1  )  =  0.0 
YK2)  =  0.0 
22   OFORMAT  (8H0OMEGA  =  ,F5.1,  5X,  11HINTERVAL  =  ,F12.10,  5X, 
115HL0AD  INERTIA  =  ,F2.1) 
WRITE  OUTPUT  TAPE  8,  22,  W(I),  TINT(I),  DL(J) 
M  =  0 

5  CALL  RKUTTA1  (2,  T,  Yl,  DELT) 
T  =  T  ♦  DELT 

FORCE  =  A*SINF{OMEGA»T) 
THEC  =  Y 1 ( 1 ) 
THECD  =  Yl (2) 
IF(T  -  DELAY)  26,  26,  30 
30    M  ■  M+  1 

IF  {M  -  100)  26,  27,  27 
27   M  =  0 

WRITE  OUTPUT  TAPE  8,  4,  T,  THEC,  THECD,  FORCE 

4  FORMAT  (F10.3,  2F20.5,  hOX.    F20.5) 
IF  (T  -  OELAY  -  FINI)  26,  19,  19 

26   IF  (BL)  5,  5,  21 

21   THECDD  =  -2-.iO«ZETA*WN»THECD  -  WN»WN*THEC  ♦  WN*WN»A»S  INF  (  W{  I )  *T  ) 

SLOPEC  =  THECDD/THECD 

SLOPEL  =  -FL/DLU) 

IF  (SLOPEC  -  SLOPEL)  6,  6,  5 

6  THEL  =  THEC 
Y3(  1)  ■  THEC 
Y2(2)  =  THECD 
Y3(2)  =  THECO 

IF  (T  -  DELAY)  16,  16,  33 
33    WRITE  OUTPUT  TAPE  8,  4,  T,  THEC,  THECD,  FORCE 
16   IF  (Y3(2))  7,  8,  8 

7  Y2( 1  )  *  (Y3(  1)  -  BL)/  RHO 
GO  TO  9 
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8  Y2(  1  )  =  Y3(  I  )/RHO 

9  CALL  RKUTTA2  (2,  T,  Y2,  DELT) 
THEM  =  Y2(  1  ) 

THEMD  =  Y2(?) 

CALL  RKUTTA3  (2,  T,  Y3,  OELT) 

THtL  =  Y5(  1) 

THELD  =  Y3(2) 

T  =  T  ♦  UELT 

J-ORGF  *  A*$INF(UMeGA*T)a 

IF  (T  -  DELAY)  *Mt  ?H,  3  I 
31    M  =  M  ♦  1 

IF  (M  -  100)  28,  29,  29 
29    M  =  0 

10  FORMAT  (F10.3,  r>F20.b) 

WRITE  OUTPUT  TAPE  8,  10,  T,  THEL,  THELD,  THEM,  THEMD,  FORCE 
IF  (T  -  DELAY  -  FIND  28,  19,  19 
28    IF  (THEL  -  RHO«IHEM)  1.5,   13,  11 

11  IF  (THEL  -  RHO»THEM  -  BL )  12,  13,  13 

12  GO  TO  9 

IF  (E)  17.  17,  13 

13  BOLVEL  =  (DM(J)  /  (DM(J)  ♦  RHO  *  RHO  ♦  DL(J)))  »(RHO  •  THEMD 
1  * (  1.0  +  E)  f  THELD  *  ((RHO  *  RHO  »  DL(J))  /DM(J)  -  E)) 

BOMVEL  =  (BOLVEL  -  E*RHO*THEMD  ♦  E*THELD)/RHO 

THELD  =  BOLVEL 

THEMD  =  BOMVEL 

Y2(2)  =  THEMD 

Y3(2)  ■  THELD 

DIFVEL    =    RHO*THEMD    -    THELD 

IF    (DIFVtL)     14,     15,     15 

14  DIFVEL  =  -DIFVEL 

15  IF  (DIFVEL  -  EPS)  17,  17,  16 

17  Yl(  1  )  =  THEL 
YK2)  =  THELD 

35   IF  ( T  -  DELAY)  34,  3U,  1 8 

18  WRITE  OUTPUT  TAPE  8,  10,  T,  THEL,  THELD,  THEM,  THEMD,  FORCE 
34    GO  TO  5 

19  CONTINUE 
END  FILE  8 
STOP 
END 

SUBROUTINE  RKUT T A  1 ( NUMBE R , XVAR , YVARS, STEP) 

DIMENSION  YVARS(30),AK(i*,30),DY(30)  ,YC(30)  ,C(4) 

COMMON  ZETA,  WN ,  A,  OMEGA,  FLt  DLOAD,  FM,  DMOTOR,  GC ,  THEL 

CM)  =  0.0 

C(2)  =  0.5 

C( 3)  =  0.5 

C(4)  =  1.0 

DO  1  1=1,4 

XC  =  XVAR  +  C(  1  )«STEP 

DO  2  J=1, NUMBER 

2  YC(J)  =YVARS(J)  *  C( I )«AK(I-1,J) 
CALL  DERIVK  XC,YC,DY) 

DO  1  J=l, NUMBER 
1  AK( I, J)  =  STtP  »  DY( J) 
1)0  3  J  =  1, NUMBER 

3  YVARS(J)  =YVARS(J)  +  ( AK(  I  , J )  +  2. *AK { 2 , J ) +2. *AK( 3 , J ) + AK ( 4 , J ) ) /6. 
RETURN 

END 

SUBROUTINE  DERIV  1  (T,  Y,  DY) 

DIMENSION  Y(2),  DY(2),  W(28) 

COMMON  ZETA,  WN,  A,  OMEGA,  FL,  DLOAD,  FM,  DMOTOR,  GC,  THEL 

DY(  1  )  =  Y(2) 

DY(2)  =  WN«WN*A»SINF(OMEGA*T)  -  WN»WN»Y(1)  -  ( 2. 0«ZETA»WN ) *Y( 2 ) 

RETURN 

END 

SUBROUTINE  RKUT TA2 (NUMBER , XVAR , YVARS, STEP) 

DIMENSION  YVARS(30),AK(4,30) ,DY( 30 ) , YC ( 30) ,C ( 4 ) 

COMMON  ZETA,  WN,  A,  OMEGA,  FL,  OLOAO,  FM,  DMOTOR,  GC,  THEL 

C(  1)=0.0 

C(2)=0.5 
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C(3)=0.S 

C(4>=1.0 

DO     1    1=1,4 

XC=XVAR    +    C( I )»STEP 

DO    2    J=l, NUMBER 

2  YC(J)    =YVARS(J)     +    C( I)»AK( 1-1, J) 
CALL    DERIV2(XC|YC,DY) 

DO.  I    4?UNUMBER 

3  YVARS(J)    =YVARS(J)    ♦     ( AK ( 1 , J )  +  2. *AK ( 2 , J ) +2. «AK ( 3, J ) ♦ AK ( 4 , J ) ) /6. 
RETURN 

END 


SUBROUTINE    DERIV    2    (I,    Y,    DY) 

DIMENSION  Y(2),  DY(2),  W(28),  DM(3) 

COMMON  ZETA,  WN,  A,  OMEGA,  FL,  OLOAD,  FM,  DMOTOR,  GC,  THEL 

DY( 1)  =  Y(2) 
0DY(2)  =  (GC  /  DMOTOR)*A»SINF(OMEGA«T)  -  (FM  /  DM0T0R)*Y(2)  - 
KGC  /  DMOTOR  )«THEL 

RETURN 

END 

SUBROUTINE  RKUTTA3 (NUMBE R , XVAR, YVARS,STEP) 

DIMENSION  YVARS(30),AK(U,30>,DY(30),YC(30>,C(4) 

COMMON  ZETA,  WN,  A,  OMEGA,  FL,  OLOAD,  FM,  DMOTOR,  GC,  THEL 

C(  1)=0.0 

C(2)=0.5 

C(3)=0.5 

C(4)=1.0 

DO  1  1=1, H 

XC=XVAR  +  C( I )»STEP 

DO  2  J= I, NUMBER 

2  YC(J)  =YVARS(J)  +  C( I)*AK(I-I,J) 
CALL  DERIV  3  { XC ,  YC,  OY ) 

DO  1  J  =  I, NUMBER 
1  AK( I, J)  =  STEP  *  DY( J) 
DO  3  J  =  1, NUMBER 

3  YVARS(J)  =YVARS(J)  ♦  ( AK (  1  , J >  +  2. *AK { 2 , J >+2. »AK ( 3, J ) +AK ( 4 , J ) ) /6. 
RETURN 

END 

SUBROUTINE  DERIV  3  (T,  Y,  OY) 
.  DIMENSION  Y(2),  DY(2),  W(28),  DL(3) 

COMMON  ZETA,  WN,  A,  OMEGA,  FL,  DLOAD,  FM,  DMOTOR,  GC,  THEL 

DY(  1  )  =  Y(2> 

DY(2)  =  -  (FL  /DLOAD)»Y(2) 

RETURN 

END 

END 

.2  .5  .8 

.8  .5  .2 

•  1  .2  .3         .4         .5         .6         .7 

.8  .9  1.0        1.2        1.4        1.7        2.0 

2.5  3.0  U.O        5.0        6.0        7.0        8.0 

9.0  10.0  20.0       30.0       40.0       50.0       60.0 

01  .005  .0033333333  .0025        .002         .0016666667 

0014285714  .00125  .0011111111  .001         .0008333333  .0007142857 

0005882353  .0005  .0004        .0003333333  .00025       .0002 

0001666667  .0001428571  .000125  .0001111111  .0001        .00005 

0000333333  .000025  .00002       .0000166667 
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APPENDIX   C 
TYPICAL  PRINT 
OF 
COMPUTER  OUTPUT 
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